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Abstract. In this paper we consider the existence of multi-sohton structures 
for the nonlinear Klein-Gordon equation I INLKGI I in M^+''. We prove that, 
independently of the unstable character of l INLKGl solitons, it is possible to 
construct a A^-soliton family of solutions to IINLKGl . of dimension 2N , globally 
well-defined in the energy space X for all large positive times. The 
method of proof involves the generalization of previous works on supercritical 
NLS and gKdV equations by Martel, Merle and the first author |3] to the 
wave case, where we replace the unstable mode associated to the linear NLKG 
operator by two generalized directions that are controlled without appealing to 
modulation theory. As a byproduct, we generalize the linear theory described 
in Grillakis-Shatah-Strauss [9] and Duyckaerts-Merle [7] to the case of boosted 
solitons, and provide new solutions to be studied using the recent Nakanishi- 
Schlag [22] theory. 



1. Introduction 

In this paper we are interested in the problem of constructing multi-soliton solutions 
for some well-known scalar field equations. Let / — f{s) be a real-valued 
function. We consider the nonlinear Klein-Gordon equation (NLKG) in d > 1, 

dttu - Au + u~ f{u) ^0, u{t,x)eR, {t,x)eRxR'^. (NLKG) 

This equation arises in Quantum Field Physics as a model for a self-interacting, 
nonlinear scalar field, invariant under Lorentz transformations (see below). 

Let F be the standard integral of /: 

F{s):^ f f{a)da. (1) 

JO 

We will assume that for some fixed constant C > 0, 

(A) If d= 1, 

[i) / is odd, and /(O) = /'(O) = 0, and 

{ii) There exists sq > such that F{sq) — -Sq > 0. 

(B) If d ^ 2, / is a pure power i/^-subcritical nonlinearity: f{u) — A|ii|?'~^M, 
where A > 0, p G (1, 3^)- 
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Prescribing / to the above class of focusing nonlinearities ensures that the corre- 
sponding Cauchy problem for (|NLKG|) is locally well-posed in H^iR"^) x H'-^iR"^), 
for any s ^ 1: we refer to Ginibre-Velo [TT] and Nakamura-Osawa [21J (when d = 2) 
for more details. 

Also under the above conditions, the Energy and Momentum (every integral is 
taken over W^) 

E[u,ut]it) = ^ J [\dtuit,x)f + \Vu{t,x)f + \u{t,x)f - F{u{t,x))] dx, (2) 
P[u,ut]{t) i J dtu{t,x)\7u{t,x) dx, (3) 
are conserved along the flow. 

Another important feature of equation (jNLKGp . still under the previous conditions, 
is the fact that it admits stationary solutions of the form u{t,x) = U{x) (i.e., with 
no dependence on t). Among them, we are interested in the ground-state Q = Q{x), 
where Q is a positive solution of the elliptic PDE 

AO-O + /(Q) = 0, Q>0, QeH^K'). (4) 

The existence of this solution goes back to Berestycki-Lions [IJ, provided the above 
conditions (in particular (ii)) hold. Additionally, it is well-known that Q is radial 
and exponentially decreasing, along with its first and second derivatives (Gidas-Ni- 
Nirenberg [8]), and unique up to definition of the origin (see Kwong [13], Serrin and 
Tang [U). 

In fact, our main result written below could be extended to more general nonlinear- 
ity under an additionnal assumption of spectral nature, namely that the linearized 
operator around Q has a standard simple spectrum. More precisely, Theorem [T] 
holds, as soon as / satisfies (i), (ii) and: 

(iii) If d = 2, |/'(s)| C|s|Pe'''*', for some p > 0, k > and all s e R. 

{iv) If d ^ 3, |/'(s)| C(l + \s\P) for some p < -^-^ and ah s e M. 

(v) —Az + z — f'{Q)z has a unique simple negative eigenvalue, and its kernel 
is given by {x ■ WQ\x e M.'^} and it is nondegenerate. 

Assumption (v) has been checked in the cases (A) and (B) (using ODE analysis), 
and is believed to hold for a wide class of functions /. (See Assumption 2]) . 

Since (jNLKGp is invariant under Lorentz boosts, we can define a boosted ground 
state (a soliton from now on) with relative velocity (5 E M''. More precisely, let 
P = (/3i, . . . ,/3d) g R"*, with 1^1 < 1 (we denote | • | the euchdian norm on R'*), the 
corresponding Lorentz boost is given by the {d + 1) x {d + 1) matrix 



As 



/ 7 -/5i7 ••• f3dl\ 

-hi 



where 7 := — , (5) 



(Z?/?-'" is the d X d rank 1 matrix with coefficient of index Pifij)- Then the 

boosted soUton with velocity (3 is 

Q^(i,x):-Q(^A^Q), (6) 

where with a shght abuse of notation we say that Q{t,x) = Qix) (namely we 
project on the last d coordinates). Also notice that (jNLKGp is invariant by space 
translation (shifts). Hence the general family of solitons is parametrized by speed 
P eR'^ and shift (translation) xq eW^: 

{t,x) ^ Qi3[t,X- Xo). 

This family is the orbit of {Q} under all the symmetries of (|NLKG| (general Lorentz 
transformation, time and space shifts), in particular it is invariant under these 
transformations: see the Appendix \K\ for further details. 

In the rest of this work, it will be convenient to work with vector data {u,dtu)'^. 
For notational purposes, we use upper-case letters to denote vector valued functions 
and lower-case letters for scalar functions (except for the scalar field Q^). 

We will work in the energy space H^iR'^) X L2(Rd) endowed with the following 
scalar product: denote U = {ui,U2)'^, V = {vi,V2)'^, we define 



{U\V) 



{Ul\vi) + {U2\V2) ^ J (uiVl + U2V2), (7) 

where {u\v) :— / uv, 



and the energy norm 

||C/||2 := {U\U) + (VuilVwi) = \\ui\\l, + \\u2\\l2. (8) 

It is well known (see e.g. Grillakis-Shatah-Strauss [9J) that ((3,0) is unstabl^^ hi 
the energy space. The instability properties of Q and solution with energy slightly 
above E[{Q, 0)] have recently been further studied by Nakanishi and Schlag, see [22| 
and subsequent works. Their ideas are further developments of the primary idea 
introduced in Duyckaerts-Merle [B], in the context of the nonlinear wave equation 
(where no Q is present). 

In this paper, we want to understand the dynamics of large, quantized energy 
solutions. More precisely, we address the question whether is it possible to construct 
a multi-soliton solution for (jNLKG|) . i.e. a solution u to (jNLKG|) defined on a semi- 
infinite interval of time, such that 

N 

u{t,x) ^''^^Qp-{t,x — Xj) as i — >■ H-oo. 

Such solutions were constructed for the nonlinear Schrodinger equation and the 
generalized Korteweg-de Vries equation, first in the L^-critical and subcritical case 
by Merle P7j, Martel [TS] and Martel-Merle [18J. These results followed from the 
stability and asymptotic stability theory that these authors developed. 

The existence of multi-solitons was then extended by Martel-Merle and the first 
author [3J to the supercritical case: in this latter case, each single soliton is 



This result was known in the physics Uterature as the Derrick's Theorem [5]. 
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unstable, hence the muUi-sohton is a highly unstable solution. It turns out that this 
is also the case for scalar field equations as (jNLKG[) . We prove that, regardless of the 
instability of the soliton, one can construct large mass multi-solitons, on the whole 
range of parameters /3i, . . . , (3^ €E Mf^ distinct, with < 1 and xi, . . . , xn G R'*- 
More precisely, the main result of this paper is the following. 

Theorem 1. Assume (A) or (B), and let /3i,/32, ■ • ■ ^Pn G M'' fee a set of different 
velocities 

Vz^j, /?, and < 1, 

and xi, X2, ■ ■ ■ , xn G shift parameters. 

Then there exist a time Tq G M, constants C > 0, and 70 > 0, only depending on 
the sets {f3j)j, {xj)j, and a solution {u,dtu) € ^([Tq, +00), x L'^{R'')) of 

^NLKG\) . globally defined for forward times and satisfying 



yt > To, 



N 



{u,dtu){t,x) - ^{Qp^,dtQi3j){t,x - Xj) 



We remark that this is the first multi-soliton result for wave-type equations. Al- 
though the nonlinear object under consideration is the same as for (NLS) for ex- 
ample, the structure of the flow is different (recall that all solitons are unstable for 
(|NLKGp , irrespective of the nonlinearity) . Hence we need to work in a more general 
framework, the one given by a matrix description of (jNLKGp . 

Let us describe the main steps of the proof. We first revisit the standard spectral 
theory of linearized operators around the soliton, and the second order derivative 
of the energy- momentum functional (see H in (fT4|) ) [9J. Since solitons are unstable 
objects, it is clear that such a theory will not be enough to describe the dynamics 
of several solitons. However, a slight variation of this functional (see in ([^5]) ) 
turns out to be the key element to study. We describe its spectrum in great detail, 
in particular we prove that this operator has three eigenvalues: the kernel zero, and 
two opposite sign eigenvalues, with associated eigenfunctions Z±. After some work 
we are able to prove a coercivity property for the operator modulo the two 
instability directions Z+ and Z_. 

The rest of the work is devoted to the study of the dynamics of small perturbations 
of the sum of N solitons, in particular how the two directions associated to Z± 
evolve. Using a topological argument, we can show the existence of suitable initial 
data for (NLKG) such that both directions remain controlled for all large positive 
time, proving the main theorem. We remark that this method is general and do not 
require the study of the linear evolution at large, but also a deep understanding 
of suitable alternative directions of the linearized operator. A nice open question 
should be the extension of this result to the nonlinear wave case, where the soliton 
decays polynomially. 

For the sake of easiness and clarity, we present the detailed computations in the one 
dimensional case d = 1. This case encompass all difficulties, the higher dimension 
case adding only indices and notational inconvenience: we will briefiy describe the 
corresponding differences at the end of each section. 
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Organization of this paper. In Section 2, we develop spectral aspects of the 
linearized flow around Qi^, which are more subtle than in the (NLS) or (gKdV) 
case. In Section 3, we present the nonlinear argument and prove Theorem [T] The 
Lyapunov functional has to be chosen carefully, as we cannot allow mixed derivatives 
of the form dtxU. For this reason, we also have to avoid the backward in time, 
compactness argument used previously in [15| 1181 13] '■ we construct the iV-soliton 
directly through a topological argument forward in time. 

Acknowlegment 
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2. Spectral theory 

In this section we describe and solve two spectral problems related to (jNLKGp . We 
will work with functions independent of time, unless specified explicitly. The main 
result of this section is Proposition [51 

2.1. Coercivity of the Hessian. First of all, we recall the structure of the Hessian 
of the energy around Q. Given Q — Q{x) ground state of Q and Qj3{x) — Q{'jx), 
where 7 = (1 — /3^)^^/^, we define the operators 

:= -a,, + Id-/'(Q), and L+ ^ -j-^d^^ +ld-f'{Q0), (9) 

so that L'p is a rescaled version of L+: 

{v (7^)) - (L+v) iix) . 

As a consequence of the Sturm-Liouville theory and the previous identity, we have 
the following spectral properties for _L+ , and therefore for . 

Lemma 1. The unbounded operator L+, defined in L^(M) with domain H'^iM), 
is self-adjoint, has a unique negative eigenvalue — Aq < (with corresponding L^- 
normalized eigenfunction ) and its kernel is spanned by dxQ- Moreover, the 
continuous spectrum is [1, +00), and is an isolated eigenvalue. 

We recall that from standard elliptic theory, is smooth, even and exponentially 
decreasing in space: there exists co > such that 

V/ceN, VxeM, 3 Cfc, |a^Q~(a;)K Cfcc-'^ol^l. (10) 

It is not difficult to check that one can take any cq satisfying < cq ^ \/l + Aq. 

Another consequence of Lemma [T] is the following fact: has a unique negative 
eigenvalue — Aq with (even) eigenfunction Qp{x) := Q~{"fx), its kernel is spanned 
by dxQji and has continuous spectrum [l,+oo). Additionally, we have 

Corollary 1. There exists t/p > such that, if v E H^{M.) satisfies (wjQ^) = 
{v\dxQi3) =0, then {L^v\v) ^ i^o||w||^i. 

We introduce now suitable matrix operators associated to the dynamics around a 
soliton. These operators will be dependent on the velocity parameter /3, but for 
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simplicity of notation, we will omit the subscript /? when there is no ambiguity. 
DefincH 

T = Tp :=-9xx+Id-/'(g^) = L+-/329^^, (11) 
(_»,;), (12, 

i - (o i°d) • <i^' 

and 

The operator H is the standard second order derivative of the functional for which 
the vector soliton R = {Q , dtQ p)'^ is an associated local minimizer. Later we 
will discuss in detail this assertion. The following Proposition describes the main 
spectral properties of H. Recall that (-I-) and (-j-) denote the symmetric bilinear 
forms on H^{R) x iy^(M) and L^(R) respectively, introduced in ([T]), and || • || is the 
energy norm defined in ([8]). 

Proposition 1. Let /3 £ R, \/3\ < 1. The matrix operator H, defined in L'^{M.) x 
L^(R) with dense domain H^(M.) x L^(IR), is a self-adjoint operator. Furthermore, 
there exist ap > 0, ^q, $_ G ^(R)^ (with exponential decay, along with their 
derivatives) such that 

i/$o=0, ($o|*-)-0, (15) 
(iJ$_|$_)<0, (16) 
and the following coercivity property holds. Let V = (fi,W2)"^ € H^{M.) x L^(R). 
Then, 

if (F|$o) {V\^-) = one has {HV\V) > ao\\Vf. (17) 

A stronger version of this result was stated by Grillakis, Shatah and Strauss in jHl 
Lemma 6.2], but the proof given there contained a gap, as noted in the errata at 
the end of [101 pa-ge 347]. As a replacement, the Proposition above (weaker than 
the original Grillakis-Shatah-Strauss result, but adequate for our purposes) was 
proposed in the errata |10| , without proof. We have not found a clear definition and 
meaning of the function in [TU], so therefore, for the convenience of the reader, 
we write the details of the proof in the following lines. 

Proof of Proposition^^ It is easy to check that H is indeed a self-adjoint operator. 
On the other hand, let V = {vi, V2Y' ■ We have from (fT4| . 



<-^i''>-(5i;'f?)|(::)) 

= (Twilwi) - [3{d^V2\vi) + f3{d^Vi\v2) + {V2\V2) 

= (L+viIwi) +/32(a,«i|a,t;i) + 2^(v2\d^vi) + {V2\V2) 
= (L^vilwi) + (/3a,wi + V2\^d^v^ + V2). (18) 
Recalling the notation of Corollary [U we define 



'Do not confuse with the transpose symbol (■)"^. 
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One can check from (HH]) that (<i>o|<I'o) 7^ and (i/$o|^o) = 0, since L^d^Qp = 0. 
Note additionally that by parity ($_|<i>o) = 0. Therefore, P3t is directly satisfied. 

We now prove ([T7|). Let V = (wi,W2)^ G H^{R) x i^(R) be satisfying the orthogo- 
nality properties 

(F|$o) = =0. 
Let us decompose wi in terms of the nonpositive spectral elements of i^, and 
I/^-orthogonally: 

vi = aQ^ + hdxQp + g, {q\Qji) = {q\dxQp) = 0. 

From the orthogonality conditions in (1171) . we have = so that a = 0, and 

hence from Corollary [1] 

{HV\V) = {L+q\q) + {pd.vi + V2\Pd^Vi + ^2) ^ i^olkllffi ^ 0. (20) 

We now argue by contradiction. Assume that there exists a normalized sequence 
= (w",W2)^ £ i7^(M) X L^(M) that satisfies the orthogonality properties 

= =0, ||V^"||2 = 1, and such that (iJV^"|V^") ^ 0. (21) 

Let us write the L-^-ortliogonal decomposition for each u": 

w" = hndxQfi + (j„, (gn|9xQ^) = 0. 

Then in view of (|20|) and ([21]) applied this time to the sequence y", g„ — > in 
and f3dxVi + ^ in i^. Now we compute 



= (T/"|$o) - J Kd^Qp - v^pdxxQp) 



= bn\\dxQl3\\L2 +13 J {hndxxQfi + dxqn)dxxQf3 + o(l) 

= bniWdxQdh + P^WdxxQdh) - I qndxxxQfi + o{l) 



Now (7„ — > in i^, so that {qn\dxxxQi3) 0, and hence 6„ as n +00. But 
in this case, v" = bndxQp + — ?► in and = PdxV^ + 0^2(1) in L^. 
Hence = + ||w2 Ili2 ^ 0, a contradiction to pT|) . 



The uniformity of ao with respect to (3 comes from the fact that 

\\dxQp\\h+l3^\\dxxQp\\h - (1 - l3^)\\dxQ\\h + - l3^?\\dxxQ\\l2 
is bounded from below by some constant c > for all |/3| < 1. 

It remains to show that < 0, namely (fT6| . If it were not the case, then 

we would have {HV\V) ^ for all V £ x . But we can check that 



Hi _ 

f3dxQi3 



'&Qp)) ^ (^^^^ 1^^) = -AollQ^lli^ < 0. 



This shows that < and concludes the proof. □ 
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2.2. Eigenfunctions of the linearized flow and Hessian. It is still unclear 
whether or not the coercivity property (jl7p - a key point in the proof of any 
stability result - is useful for us, since solitons are actually unstable. It turns out 
that for our purposes, we need a different version of Proposition[Tl for the linearized 
operator of the flow around Q. In order to state such a result, we introduce some 
additional notation. 

Let /? G K, |/3| < 1 be a Lorentz parameter, and consider the operators T , J, L and 
H defined in Let 

^ = ^W) = JL= (^^ , (22) 



1!^ L ]=-HJ- (23) 



and 

— I 

Id -I3d,^ 

Concerning this last operator, we prove the following result. 

Lemma 2. Let ;3 e K, < 1, 7 = (1 - (3'^)-^/'^ and Aq from LemmaUl Th ere are 
functions Zq — Z^ p, and Z± — Z±fj, with components exponentially decreasing in 
space, satisfying the spectral equations 

.J^Z^^O, and ^Z± = ±^^Z±. (24) 

7 

Moreover, by the nondegeneracy of the kernel spanned by $0, we can assume <f>o = 

Proof. The proof is similar to that of [S] . In particular, we obtain explicit expressions 
for Zo and Z± in the following lines. 

The eigenvalue problem ,^Z = \Z reads now, with Z{x) = (^1(72;), ^2(72;))"^", 

TZ2 + iPZ[ + \Zi =0, Zi - ^7Z^ - \Z2 = 0, (25) 
Replacing Zi in the first equation above, we get (recall that Q[j{x) = Q{'yx)), 

-^^z'i + Z2- nQp)Z2 + hi^z'2 + hz'i) + Khz'2 + \Z2) = 0, 

namely 

- Z!{ + Z2- f{Qp)Z2 + 2I3^\Z'2 - -\^Z2. (26) 
Performing the transformation Z2{s) :— Z2{s)e^'*'^'' , where s G M, we get 

-Z!^ + Z2- f'{Q)Z2 = -(/327' + 1)A'^2 = -\^l''Z2. 

Therefore, by virtue of Lemma[T]we can take Z2 = Q^{s) and A±7 = ±v^, where 
— Ao < is the first eigenvalue of the standard Schrodinger operator L^, defined in 
dH). Thus, 

^±,2(s)=Q-(s)e±^^^ 
Note that from pTI)l . Z±^2 decreases exponentially at both sides of the origin, since 
< 1 and PyfX^ - Vl + Ao < . 

From ([25| . we have 

Z±As) = fiiZ'^^{s) + \±Z±a{s) 



By the same reasons as above, Z±^i is an exponentially decreasing function. From 
these identities, we have 



V Q~iix) J 

V Qfi (^) / 

Now, we consider the computation of Zq. Replacing A = in (|26p . we can choose 
ZoAs) = Q'is), and Zo,i = PlQ"{s), 

from which we get 

It is clear that ,^Zq — 0. Similarly, we have M'Z± — ±-^^-^Z±, which proves 
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(HI. □ 

In order to prove Proposition [21 we need to prove the existence of two additional 
functions, both associated to Z±. 

Lemma 3. There exist unique functions Y±, with components exponentially de- 
creasing in space, such that 

HY± = Z±, ($o|i^±) = 0. 
Moreover, Y± satisfy the additional orthogonality conditions (Y±\HY±) = 0. 

Proof. Let us prove the existence of Y±. It is well-known that a necessary and 
sufficient condition for existence is the following condition: it suffices to check that 
Z± are orthogonal to (f>Oj the generator of the kernel of H . Indeed, we have from 
([24| . ([23| . the self-adjointedness of H and Proposition [U 

($o|^±) = ±^{W^Z±) = T^($o|i? J^±) = 0. 
V Ao V Aq 

However, we need some additional estimates on Y±. In what follows, we write 

down explicitly the equation HY± — Z±. It is not difficult to check that Y± = 

(y±_i, y±_2)"^ satisfies the equations 

TY±s - /3(r±,2). = ^±,1, /3(F±,i). + >±.2 = ^±,2. 
Replacing the second equation in the first one, we get (cf. @) 

=/3(^±,2)x+^±,l. 

Note that {l3{Z±^2)x + Z±^i\dxQp) = 0. Therefore, Y±^i exists and it is exponentially 
decreasing, with the same rate as Z±^i and Z±_2- A similar conclusion follows for 
Y±,2- 

Since Y± is unique modulo the addition of a constant times $0, it is clear that we 
can choose Y± such that (<I'o|^±) = 0. On the other hand, from Lemma [21 

{Y±\HY±) - (r±|Z±> - ±^(y±|^^Z±) = T^{HY±\JZ±) 

V AO V Ao 
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= T^{Z±\JZ±)^Q. 
V Ao 

□ 

The main result of this section is the following alternative to Proposition [TJ 

Proposition 2. There exists fin > such that the following holds. Let V E x 
such that ($o|y) = 0. Then 

{HV\V) > ^ioWVf --\{Z+\V)^ + {Z^VY 
/io L 

Proof . It is enough to prove that (^olV") = {Z+\V) = = imply 

{HV\V) > ^L^WVf , 

for some \xq > 0, independently of V. In order to prove this assertion, we first 
assume /? 7^ and decompose orthogonally V and Y± (cf. the previous Lemma) as 
follows 

F = y + a_$_, r± = y± + (5o$o + ^±*-, (29) 

(recaU {V\^q) = 0) with 

(y|$o> = (f±|$o> = - (i±|i/$-> - 0. 

Since ($o|$-) = ($o|^) = ($o|>±> = and ($_|iJ<I>_) < 0, it is clear that 
c^a = = and a_,(5± are well-defined. Moreover, 

Claim. For all /3 e (— 1, 1)\{0}, 1+ and Y- are linearly independent as L^(M)^ 
vector-valued functions with real coefficients. 

Indeed, to see this, assume that there is A 7^ such that y+ = XY-. Then, from 
the previous decomposition and Lemma [3l 

Z+ - AZ_ = H{Y+ - Ay_) = {S+ - A(5_)iJ$_. (30) 

This identity contradicts ^27} and p^ . which establish that Z+ and Z_ have 
essentially different rates of decay at infinity, different to that of for all 

/3 7^ 0, which makes ([5U)) impossible. 

The analysis is now similar to that in [7, Lemma 5.2]. We have from (|29p . 

{HV\V) ^ + + = {HV\V) + (31) 

On the other hand, since {Z±\V) — 0, we have from Lemma [31 

= {Y±\HV) = {Y± + S±<^>_\HV + a_H<i>_) = (f±|iJV^) + a_5±(i7$_|$_). 
Similarly, 

= (if y± I y± ) = (if f ± I f± ) + 4 (if $_ I ) . 

We get then 

(iiT/lF) ^ (ii^lF) - . (32) 



Consider 



sup 

WeSpa.n{Y+S-)\{0} 



'{HY+\Y+){HY^\Y^) 
(y+|iiM^) (f_|iiT4^) 



'HY+\Y+){HW\W) ^^/{HY^\Y^){HW\W) 
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Recall {H ■ \ -) is positive definite on Span(<i>o, $-)^- Hence apply Cauchy Schwarz's 
inequality to both terms of the product: it transpires that a ^ 1. Furthermore, if 
a — 1 (as Span(y+,F_) is finite dimensional), there exists W of norm 1 such that 
both terms are in the equality case in the Cauchy-Schwarz inequality, i.e W and 
are linearly dependent, and W and Y^ are also linearly dependent. But it would 
then follow that 1+ and Y- are linearly dependent, a contradiction the above claim. 
This proves a < 1. 

Now using i?-orthogonal decomposition on Span($o, ^-)^, we deduce that 



\/W& Span($o,*-)^ 
By (131D and (HIl), we get 



{Y4HW){Y+\HW) 



{HY+\Y+){HY^\Y^ 



a{HW\W). 



{HV\V) ^ (1 - a){HV\V) ^ ao{l - a)\\V\\'^ > 0. 

and so dSI]) implies {HV\V) ^ 
We then conclude that, for C ~ tt^t- max ( — 



C{HV\V) ^ C(l - a){HV\V) 

^ 2||t/|p + 2a2_||$_||2 ^ \\V + a-^-\\^ = \\V\\^. □ 
Finally, if /? = 0, we proceed as follows. First of all, we have from (P7| and 

so that {^o\V) ^ {Z±\V) = imply (vilQ') = (vilQ-) = (1^2 = 0, where 
V — {vi,V2)'^- Additionally, from ^TE\\ . we have for some e > small and fixed, 

{HV\V) - {L+vi\vi) + Wd^vi + V2Wd^vi + V2) 

= iL+V,\v,) + (3^d,V,\\^ + 2P{d,V,\v2) + \\V2\\^ 

> ^^olkillffi - Y^\\d,vi\\^ + e\\v2\\^ 

2.3. Extension to higher dimensions. The equivalent of Lemma [T] (and there- 
fore assumption (iv) of the Introduction) in dimension d > 2 has the form 

Lemma 4. Assume d ^ 2 and assumption (B) holds. has exactly one negative 

eigenvalue, and its kernel is spanned by {dxiQ)i=i,....d- Its continuous spectrum is 
[l,+oo). 

Proof. See Maris [M] and McLeod [l6]. □ 

As mentioned in the Introduction, this result is open for general nonlinearity /. In 
that case, we need to assume that it holds, i.e. assumption (v). 
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The null directions for H are now the d-dimensional vector space spanned by the 
functions $o,i = ( r^) ■ the proof of Lemma [21 one should rather per- 



form the transformation Z2 = Z2e^^^^'^ . The rest of the arguments is dimension 
insensitive. 



3. Proof of the Main Theorem 

In this section we prove Theorem [TJ Again, we will give a detailed proof in the one 
dimensional case d = 1, and point out how to extend the proof in higher dimension, 
which is done in a similar fashion as in |15j . 

3.1. The topological argument. We continue with the same notation as in the 
previous section. In particular, we fix /? G (—1,1) and consider now the time- 
dependent, boosted soliton given by 

Qp{t, x) = Q{^{x ~ pt)), 7 = (1 - p^r^'^- 

Additionally, we suppose given iV different velocities . . . ,f3M G (—1, 1), already 
arranged in such a way that 

- 1< /3i < /32 < . . . < /3iv < 1, (33) 

and N translation parameters xi, . . . ,xm € R, such that Qp^{t,x — Xj) is the 
associated soliton solution of velocity /3j and shift Xj, j = I, . . . , N. 

Finally, we introduce some notation. Given ^ a real Banach space, x £ B and 
r ^ 0, we denote 

Bfij{x,r) ^{yeJ^l \\x~y\\fij s: r} 

the closed ball in centered at x of radius r and || • ||^ is the associated Banach 
norm on 

Lemma 5 (Modulation). There exist Lq > and Eq > such that the following 
holds. For any L ^ Lq and < e < eo> io G ifU€ iJ^(K) x L^(R) is sufficiently 
near a sum of solitons whose centers are sufficiently far apart, 



min {\yj - y,\\ i ^ j} ^ L, 



then there exist shifts yj — yj{f3j,to) such that if we define 



then for some C > 



i?,(x):= [d%yO,x~y,), (34) 

N 

R{x):=Y,R,{x), (35) 

Vix) ■.= Uix) - R{x), (36) 



\\V\\^Ce, and {V\{Rj)^) = 0. (37) 
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Also, the map U ^->■ (V, {yj)j) is a'^^-dijfeomorphism around^^^^ (yd^Q ) i^Oi^— 

Uj). In such case, we say that U can be modulated into (V, {yj)j)- 

Proof. This is the classical modulation result, stated as in [3, Lemma 2]. See [26l[27] 
for more details. □ 

In what follows, we introduce additional notation. We assume that U can be mod- 
ulated into {V, {yj)j)- 

For any j ~ 1, . . . , N (cf. Proposition [1] and Lemma [2] for the definitions), let 



(38) 



(39) 



I ^±j(s) := 2'±(7j(s - y^)), Zoj{s) := Zo{-fj{s - y,)), 
\$Oj(s) :=$o(7j(s-y,)), $-j(s) :=$-(7j(s-yj)), 
where 7j :— (1 — /3|)^-^/^, and 

|a±,, (F|Z±,,), 

along with the vectors 

a+ = (a+j)j, a- = (a_j)j, = (ao,j)j, and y ^ (40) 
Finally, we fix a constant 70 given by 

70 := imin{7i,72,...,7Ar}min{/3i,/32 -/3i,...,^Ar - Pn-i} > 0. (41) 

Assume now e G (0, £0) and L > Lq, where eg and Lq are obtained in by Lemma [S] 
Given < G M, let us consider the centers 

where the velocities Pj and the shifts xj are given by It is clear that there 

exists To G Hi such that, for all t > Tq, the yj satisfy 

minjlyj -yt\\i ^ j} > L. 

From now on. we fix i > Tq. Consider the corresponding sum of solitons R(t,x) 
associated to these parameters, namely 

R{t,x) f2 ^^(*' ^) = E ( d%. ) " y^^- 

j=l j=l \ 

Then, according to Lemma[5l if C/ S H^(R) x L'^{R) satisfies \\U-R(t)\\ ^ e, then U 
can be modulated. Moreover, up to increasing Tq, we can assume that e"'^""^" < Eq. 
Thus we can define our shrinking set. 

Definition 1 (Shrinking set y{t)). For t ^ Tq, we define the set 

^(i)cSffixL2(i?(t),eo) 
in the following way: U E 'fit) if and only if U can be modulated into (V, y) where 
(cf. dMl) and (gni)) 

AT 

v{t) = u-Y,m), 

with 

\\V{t)\\ < e-^«*, \y^{t)-p,t-x,\ ^ e-^«*, (43) 

13 



Ila+ll,. <e-3^«*/2, ||a_||p ^;e-3^«*/^ ||ao||p ^ 6"'^°*^'. (44) 

Definition 2. We denote hy if = {u, dfu)^ the flow of the (NLKG) equation, that 
is, given toGR and Uq G i/i(M) x ^^(R), 

th^(/7(to,i,C/o) (45) 
is the solution to (NLKG) with initial data Uq at time to (with values in x L^). 

Consider an initial data Uq at time Tq such that Uq E YiTo). Due to the blow- 
up criterion for (NLKG), and the fact that R{t) defined in (|42l) is bounded in 
i/^(R) X L^(M), we have that (p{Tq, t, Uq) is defined at least as long as it belongs to 
BH^xL^{R{t)^ !)■ Ill particular, ip{T(),t, Uq) does not blow-up as long as it belongs 
to 'V{t), and we can define the exit time 

T*(C/o) :=sup{T^To |VtG[To,T], ^(Tq, t, C/q) G r(t)} . 

Notice that we can have T*{Uo) = Tq. Our goal is to find Uq € l^(To) such that 
T*{Uo) - +00. 

In order to show such an assertion, we will only consider some very specific initial 
data, namely Uq G fiTo) such that (see pO|) ) 

• a_(To) = ao(ro), and 

• a+{To) G BR«(0,e-3-ro7'o/2). 

These conditions are satisfied due to the almost orthogonality of Z±j, Zqj, and is 
the content of the following 

Lemma 6 (Modulated initial data). Let Tq > Q he large enough. Given a+ = 
{aj^j)j G i?RN(0,l), there exist a unique rescaled Uq := 0(a+) G YiTo) such that 
Uq can be modulated into (Vo,y) and the associated parameters (|40p satisfy 



a+{To) - e-^-'°^°/^a+, a_(To) = ao(ro) = 0. (46) 

Moreover, 

\\Va\\ <Ce-^^°^°/^. (47) 

Proof. As in [Si Lemma 3]. □ 



If T*{Uo) is finite, by maximality, we also have that for the function ip(TQ,T*, Uq), 
at least one of the inequalities in the definition of y{T*) is actually an equality. 
It turns out that the equality is achieved by a+{To) only, and that the rescaled 
quantity e'^'^°'^ /^a+(t) is transverse to the sphere at t = T*. This is at the heart 
of the proof and is the content of the following 

Proposition 3. Let a+ G i?jjiv(0, 1), and assume that its corresponding maximal 
exit time is finite: 

T* = r*(e(a+)) < -foo. 
Recall, for all t G [To,T*], the associated modulation {V(t),y(t)) o/(p(t, Tq, 0(a+)), 
defined in Then, for all t G [To,T*], 

\\Vm ^ \y,{t) - p,t-x,\ ^ ie-^°*, (48) 

and 

|la+(T*)||,. =e-3''o^-/2. (50) 
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Furthermore, a^(T*) is transverse to the sphere, i.e., 

d 



> 0. 

For the sake of continuity, we postpone the proof of Proposition [3] until the next 
section, and conclude the proof of Theorem [1] here, assuming Proposition [3l 

Let us state a few direct consequences of Proposition [3l (their proofs will also be 
done in the next section [ 



Corollary 2. We have the following properties. 

(1) The set of initial data which give rise to solutions which exit in finite time 

n := {a+ e Bk«(0, 1) I T*(e(a+)) < +c»} 
is open {in B^n{0, 1)). 

(2) The map 9 a+ i— > r*(0(o+)) G M is continuous. 

(3) The exit is instantaneous on the sphere: 

if ||a+||f2=:l, then T*{a+)^Ta. (51) 

We are now in a position to complete the proof of Theorem [1] 

End of the proof of Theorem QJ We argue by contradiction. Assume that all possi- 
ble a+ e i3Riv(0, 1) give rise to initial data Uq — 6(a+) G yiTi^) and corresponding 
solutions Lp{tQ,t,U()) that exit '^{t) in finite time, i.e. 

assume that il = B^n{{), 1). (52) 

Given Uq G YITq), we denote $([/o) the rescaled projection of the exit spot 

$([/o)=e3^°^*(^°)/2a+(r*([/o)), 

so that $(C/o) G Sr"(0, 1). Let us finally consider the rescaled projection of the 
exit spot ^ , defined as follows: 

: Br«(0, 1) ^ BKiv(0, 1), a+ ^ *(a+) = $ o e(o+). 

Corollary [2] then translates into the following properties for VP: 

• ^ : -Br]v(0, 1) — >■ B^n{(), 1) is continuous (as T* and $); 

• If ||a+||£2 = 1, *(a+) = a+ (cf. 1^ and dM])). 

These two affirmations contradict the Brouwer's Theorem. Hence our assumption 
([5^ is wrong, and there exists a+ such that the solution U{t) = tf{TQ,t,<d{a+)) 
satisfies U{t) G '^(t) for all t ^ Tg. In particular, U{t) satisfies the conditions of 
Theorem [H □ 

3.2. Bootstrap estimates. This Section is devoted to the last remaining results: 
Proposition [3] and Corollary [2] 

Proof of Proposition\^ Step 1. First, we introduce some notation. Consider the 
fiow ip{t) = (p(t,ro, e(a+)) given by Proposition [3l and valid for all t G [To,r*]. 
From Lemma [SI we have 

^{t) = Kit) + V{t), (53) 

where 

N 

i?(t,x) =^i?,(i,x), Rj{t,x) = {Qp^,dtQpf{x-y,{t)), (54) 
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(55) 



and 



Additionally, from the equation satisfied by Lp, we have 



Id 

52 -Id 





/(") 



where tp = iu^utf" . Replacing the decomposition (j53p . we have 



Vt 



Id 



with = ^{fi.j) defined in ([22]), 



i?e(i) 



Id 

52 -Id 



R — Rt 



V + Re{t) :^ + Re{t), 





f{u)-nQpM 



(56) 



-|:4(t)a.(_,X)J 



First of all, note that from (1191) we have 



XX 

is a constant C > 0, independent of time, such that 



$o,fc- Additionally, there 



N 



N 



k=l 



fe=i 

fc=i fe=i 

N N 



fe=l 



fc=l 



-/'(E 

TV 



^ C(e-37ot^_^2), (57) 

(see (HI])). Finally, if we take the scalar product of ([5S)) with {Rj)x, then the or- 
thogonality (I57|) (coming from modulation) leads to the estimate 

|i;(i)KC(rWII+e-^^°*), (58) 

valid for all j = 1, . . . , TV. 

Step 2. Control of degenerate directions. The next step of the proof is to con- 
sider the dynamics of the associated scalar products a±j{t) and ao,j{t) introduced 
in 

Lemma 7. Let a±j{t) and aoj{t) be as defined in f39\}. There is a constant C > 0, 
independent of To and T* , such that for all t G [To, T*], 



Ao 



a'±,,{t)T—a±,,{t) 
I3 



^ C\\V{t)f + Ce-^^°\ 



(59) 
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and 

Wo.jit)] ^ C\\V{t)\\^ + Ce-^^°*. (60) 

Proof. We prove the case of a^j{t). The other cases are similar. We compute the 
time derivative of a+j- using ([Ml) and ([55)) . and we choose 70 > as smah as needed, 
but fixed. 



= -y;-(t)((^+,,).|y(t)> + (Z+Mit)) 



N 



k=l 

+ 0(||y(<)f + e-37ot). 

From Lemma[3]we have (<i>oj|Z+j) — 0. Therefore, since ^* — /3jdx — J^j, where 
J'fj := .y/f{Qp^) (cf. ([231)), we have from Lemma[2]and ((58)) . 



= ^^+At) + oiw^mm + \\v{t)r + e-^^°*) 



^%+,,(<) + 0(||F(t)f + e-37ot). □ 



Ij 

Step 3. Lyapunov functional. Let Lq > be a large constant to be chosen later. 
Let (0j)j=i,...,Ar be a partition of the unity of M placed at the midpoint between 
two solitons. More precisely, let cj) S C°°(]R), (j>' > and lim_oo <P — 0, lim+oo 0=1. 
We hav^, for all L > Lq, 

^,(^,a:) = l, </,^.(t,a:) ^^j, (61) 

where nij :— \{l3j + /3j-i), with j = 2, . . . , iV — 1, and mi := —00, m^r = +00. We 
introduce the j-th portion of momentum 

PjYp]{t) := ^ / <t>]UtUx dx, ip = {u,utf , (62) 



2 . 

and the modified Lyapunov functional 

N 

^M(t) := E[ip]{t) + 2^/3,P,-M(t), (63) 

i=i 

with being the energy defined in ([2]). Our first result is a suitable decomposition 
of J?[m] around the multi-soliton solution. 

Lemma 8. Let V{t) = {vi[t),V2{t)Y be the error function defined in Proposition 
\3i There is a positive constant C > such that 

N 

^M(i) - Y.(HjV\V) < C||y||^ + -e-2^°*, (64) 

where 



■^Do not confuse the constant L in I I13I I with the operator L in 11131 . 
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Proof. From the decomposition 

ifiit) = {u,Ut){t) = {Ri,R2f + {Vi,V2f{t), 

we have 



(66) 



•I J=l 

N 

+ / [R2V2 + {Rl)xivi)a: + RlVi - f{Rl)Vi + ^ djiMvi)^: + V2{Rl)x)(pj] 

- J {F{R, + vi) - F{R,) - f{R,)v, - ^f{Ri)vl) 
=: Ii+l2+h + h- 

Let us consider the term h . Since R2 = - J^f^i ft {Qffj )x and (.Ri )x = Z^^i (<3/3j )x , 
one has 

1 ^ r 



N 



Now we consider /2. Integrating by parts, we have 



= / ^;2[-R2 + (^l)xEft <^: 



AT 



JV 



- / t^l[(^l)xx-^l + /(^l) + (^2)xEft-'^j] -Eft/ ^1^2(<^i)x. 

Note that 

N N N N N 

R2 + (^1). Eft-^i = J2{-MQ0>.)- + w^Jx Eft-^^i = E(Q/3J- Eft'^j-- 

J=l fe=l j=l fc=l JT^fe 

Hence 



/ 



V2 



N 



R2 + iRi).^PAj\ =0(e-3^°*). 
On the other hand, Therefore, 



/ 



N 



Vi 



(El),, - i?i + /(i?i) + (i?2)x Eft<^jJ = 0(6-3-^°*). 
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< C'||i;i||i2(R)p/ 



Finally, 

' ^ r 

Gathering the above estimates, we get 

1/2 1 s; Ce-3^«*. 

Let us consider the integral . Since = 1 , we have 

N 



270* 



N 

Therefore, 

h = \{H,V\V)+0{e''-^'^'). 
Finally, we consider 14. It is not difhcult to check that 

Collecting the above results, we get finally ([M)) . □ 
Our next result describes the variation of the momentum Pj. 

Lemma 9. There exists C > independent 0/ time and L, such that for all t € 
[To,T*], 

|P,M(t)-P,M(To)K|e-2'^«*. (67) 
Proof. A simple computation using (INLKGP shows that 

^ ^ F{u){^j), + ]- j utu,{(j)j)t. (68) 



2 ./ ' '^-^-i'-^ 2 

Indeed, one has 

dtPj[ip\{t) ^ \ I UtUx{4>j)t + \ ( UtUtx4>j + \ I UuUx4'j 



2 

^ / UtUx{(t>j)t + ^ I {ul)x(t)j + ^ / (""l^^r -U + f{u))u^ 



= ^ J utUx{(j)j)t Uti(j}j)a; - ^ J{ul-u^ + 2F{u)){(j)j)x, 

as desired. Now, from the decomposition (j66|) we replace above to obtain (compare 
with ([57])) 



g-37ot 

^ ( —r- + / "2(0.)x + / {vifA4>j). + / vi{cp,)x + / F(i;i)(0,). 
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From the smallness condition of v, we get finally 

C 
L 

as desired. The conclusion follows after integration in time. □ 
The previous Lemma and the energy conservation law imply the following 
Corollary 3. There exists C > independent of time and L > such that, for all 

|^M(t)-^M(To)K^e-27ot. (69) 

Now we use the coercivity associated to Hj . A standard localization argument (see 
e.g. [2n])j Proposition [2] and (|39p give 

^ 1 

for an independent constant vq > 0. From this coercivity estimate, using ([M)) and 
([M)) . the initial bound (|T7| . and bounding the terms in a± by ([33]), we get that for 
some C > 

Vie [Tq,T*1 \\v{t)\\ -^e~^°* + Ce-3/27ot_ 

vL 

Therefore, for L ^ AC^, we improve the first condition in P5)) . to get We can 
now integrate of the modulation equation ((58)) for x'^ {t) we get the second estimates 
in (|48p (by increasing L is necessary). 

Now, using (|59 p -()60 p and integrating in time, we improve in a similar way the 
conditions in (|331), to obtain (H^ . 
In conclusion, ([5n|) must be satisfied. 

Step 4. Transversality. For notation, let ^(a+,i) := e^'''"*||a+(t)|l|2 . Using the 
expansion (|59p . we compute 

^(a+,i) = ^e3^°*a;^^.(0a+j(i) +37o^(a+,0 



N 



= e3^°*yA^f]^±A+O(e370*(|lF(t)|li.+e-2^°*)ll«+Wll.0 
+ 37o^(a+,t) 

> (2co + 37o)^(a+,<) + Oie^''»\\\Vmh + e-2^°*)||a+(<)||,0, 
where cq = iVAj^min,{l/7,} > 0. Now, at time T* = T*(a+), \\V{T*)\\l2 
0(e-'^o^'), whereas ||a+(T*)||^2 = e^37oTV2^ j g. ^(a+,r*) = 1 , hence 

d 



^ (2co + 37o)+0(e-T«^'/2) 

t=T*(a+) 

Choosing Tq larger if necessary, and as T* ^ Tq for all o+, we get 



|^(a„.) 



Co + 370 > 0. (70) 
This concludes the proof of Proposition |3l □ 
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We end this paragraph with the proof of Corollary 

Proof of Corollary [H Let us now show that is open and that the mapping a+ i— >■ 
T*(a+) is continuous. Let a+ e n. We recall that ^{a+,t) = e^'^^'^Wa+it^^. By 
(ffO)) . for all e > small, there exists 6 > such that 

• ■yV{a+,T*{a+) +e)> 1 + 6, and 

• for aU t £ [To,T*{a+) - e] (possibly empty), J^{a+,t) <1~S. 

By continuity of the flow of the (|NLKG[) equation, it follows that there exists 77 > 
such that the following holds. For all a+ G i3RN(0,l) such that ||o+ — o+|| ^ r/, 
then \^{a+,t)-jK{a+,t)\ ^ S/2 for allte [To, T*(a+) +e]. In particular, a+ £ n 
and 

r*(a+) - e ^ T*{a+) s: T*(a+) + e. 
This exactly means that ft contains a neighbourhood of a+ , hence is open, and that 
a+ 1—^ r*(a+) continuous. 

Finally, let us show that the exit is instantaneous on the sphere. If ||a+||£2 = 1, 
then ^(o+,To) = 1, hence by ((701), ^(a+,t) > 1 for aU t > Tq in a neighborhood 
of Tq. This means that r*(a+) = Tq. □ 

3.3. Extension to higher dimension. The main part of the proof remains un- 
changed. One has to work only for the definition of the Lyapunov functional. The 
key point is to notice that one can find a suitable direction as in |15) . The set 

^ = {/3 G M'' I Vj, /3./3j=0}, 

is of zero measure: let ^ ^ up to rescaling, we can assume \P\ = 1. Without loss 
of generality we can assume that the indexes j satisfy 

-1< • /3i) < • /32) < • • • < • M < 1. 

We use again the Id cut-off function cj) defined at Step 3 of the previous to define 
the new cut-off functions 

Then all the computations of Step 3 of Section 13.21 follow unchanged. We refer to 
[15] (Claim 1 and what follows) for further details. 

Appendix A. The orbit of Q under general Lorentz transformations 

In this appendix we prove that the orbit of Q under the group generated by space 
and time translations, and general Lorentz transforms is 

^ := {{t, x) ^ Qp{t, x-xo)\ /?, xo G < 1}. 

We recall that we consider Q as a function of time with the slight abuse of notation 
Q{t,x) = Q{x). 

The map /3 A^ (see ([S])) is a group homomorphism from (i?jjd (0, 1), ©) to 
(A/c;+i(M), o), where © denotes Einstein's velocity addition 

1 / x ■ y r / X ■ y W 

x®y^ — U+ -ri2"2/+ Vl - \v\ \x- 

l + x-y \ \y\^ \ \yY J j 

In particular, A-^jA^ — Id^+i. 
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A general Lorentz transform is an element of 0{l,d) 
written in the form 

/I • • • 0\ 




x: 0{d), hence can be 



A 



U.fi 



U 



Vo 



A^, where U G SO{d), i.e UU' 



As Q{x) is radially symmetric, it follows that 



Q A 



Qp{x) 



hence the orbit of {Q} under general Lorentz transform is simply {Qp]^ € R*^}. We 
now want to parametrize the other invariances of (|NLKGp . that is time and space 
shifts. Fortunately, the former reduce to the latter. Indeed, notice that 











V 



T — 1 

so that Idd +4:^/3/3^ 



(l 




Idd-i 



Vo 



(Here ^ indicates similarity of matrices). In particular, Id^ is invertible. 



Then time translations for Qp can be rethought as an adequate space shift: 



t + to 



= g A 



7 



g 



to 



It follows that is stable through all general Lorentz transform, time and space 
shifts, hence it is the orbit of Q through the group generated by these transforma- 
tions. 
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